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Abstract 

Recently, we developed a method for calculating the lifetime of a 
particle inside a magnetic trap with respect to spin flips, as a first step 
in our efforts to understand the quantum-mechanics of magnetic traps. 
The ID toy model that was used in this study was physically unrealis- 
tic because the magnetic field was not curl-free. Here, we study, both 
classically and quantum-mechanically, the problem of a neutral par- 
ticle with spin S, mass m and magnetic moment /i, moving in 3D in 
an inhomogeneous magnetic field corresponding to traps of the loff^e- 
Pritchard, 'clover-leaf and 'baseball' type. Defining by ujp, iVz and 
iOr the precessional, the axial and the lateral vibrational frequencies, 
respectively, of the particle in the adiabatic potential Ve/f = /^|B|, 
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we find classically the region in the {uJr/cOp)-{uJz/^p) plane where the 
particle is trapped. 

Quantum-mechanically, we study the problem of a spin-one par- 
ticle in the same field. Treating oj^/oop and ooz/ojp as small parame- 
ters for the perturbation from the adiabatic Hamiltonian, we derive 
a closed- form expression for the transition rate l/Tesc of the particle 
from its trapped ground-state. In the extreme cases the expression for 
l/^esc reduces to 



Ai^LOr exp 



Sy2'K yjujpUJi exp 



2ujp 

UJr 
2uj^ 

UJi 



; for ujp » UJr » uJz 



; for ujp ^ UJr 



UJz 



UJi 



UJr 



(uJzV" 


2ujp 


— exp 






UJz _ 



; for cjp » » UJr 



1 Introduction. 

1.1 Magnetic traps for neutral particles. 

Recently there has been rapid progress in techniques for trapping samples 
of neutral atoms at elevated densities and extremely low temperatures. The 
development of magnetic and optical traps for atoms has proceeded in parallel 
in recent years, in order to attain higher densities and lower temperatures 

0, |, ^, K . We should note here that traps for neutral particles have been 
around much longer than their realizations for neutral atoms might suggest, 
and the seminal papers for neutral particles trapping as applied to neutrons 
and plasmas date from the sixties and seventies. Many of these papers are 
referenced by the authors of Refs.jjl], @, In this paper we concentrate 
on the study of magnetic traps. Such traps exploit the interaction of the 
magnetic moment of the atom with the inhomogeneous magnetic field to 
provide spatial confinement. 

Microscopic particles are not the only candidates for magnetic traps. In 
fact, a vivid demonstration of trapping large scale objects is the hovering 
magnetic top [^, |^, ^. This ingenious magnetic device, which hovers in 
mid-air for about 2 minutes, has been studied in the past few years by several 



authors |1C, 11, 12, 13, 14, 15 



2 



1.2 Qualitative description. 

The physical mechanism underlying the operation of magnetic traps is the 
adiabatic principle. The common way to describe their operation is in terms 
of classical mechanics: As the particle is released into the trap, its magnetic 
moment points antiparallel to the direction of the magnetic field. Inside 
the trap, the particle experiences translation oscillations with vibrational 
frequencies u^ib which are small compared to its precession frequency ujprec- 
Under this condition the spin of the particle may be considered as experi- 
encing a slowly rotating magnetic field. Thus, the spin precesses around the 
local direction of the magnetic field B (adiabatic approximation) and, on 
the average, its magnetic moment fi points antiparallel to the local magnetic 
field lines. Hence, the magnetic energy, which is normally given by — /i ■ B, 
is now given (for small precession angle) by |B|. Thus, the overall effective 
potential seen by the particle is 

Veff^fi\B\. (1) 

In the adiabatic approximation, the spin degree of freedom is rigidly coupled 
to the translational degrees of freedom, and is already incorporated in Eq. (|ip 
such that the particle may be considered as having only translational degrees 
of freedom. When the strength of the magnetic field possesses a minimum, 
the effective potential becomes attractive near that minimum, and the whole 
apparatus acts as a trap. 

As mentioned above, the adiabatic approximation holds whenever utprec 3> 
Uyib. As Wprec IS inversely proportional to the spin, this inequality can be 
satisfied provided that the spin of the particle is sufficiently small. If, on 
the other hand, the spin of the particle is too large, it cannot respond fast 
enough to the changes of the direction of the magnetic field. In this limit 
i^prec ^ ^vib, the spiu has to be considered as fixed in space and, according 
to Earnshaw's theorem ||T6|, becomes unstable against translations. Note 
also that Uprec is proportional to the field |B|. To prevent Uprec of becoming 
too small, resulting in spin-flips (Majorana transitions), most magnetostatic 
traps include a bias field, so that the effective potential Veff possesses a 
nonvanishing minimum. 
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1.3 The purpose and structure of this paper. 

The discussion of magnetic traps in the hterature is, almost entirely, done 
in terms of classical mechanics. In microscopic systems, however, quantum 
effects become important, giving rise to a finite lifetime of the particle within 
the trap. This requires a quantum-mechanical treatment [0. An even more 
interesting issue is the understanding of how the classical and the quantum 
descriptions of a given system are related. It is important to note here that 
there are two mechanisms by which the particle can escape from the trap: 
The first one is the usual tunneling of the particle from the trap, without 
a change of its spin state, to regions where the magnetic field decreases to 
zero. The time scale for this process can be evaluated by standard methods. 
The second way the particle can escape from the trap is by flipping its spin 
state (Majorana transitions). This process, which is different from the first 
one because there is no potential barrier, is the subject of this paper. 

As a first step in our efforts to understand the quantum-mechanics of 
magnetic traps, we recently developed a method for calculating the lifetime 
of a particle inside a magnetic trap with respect to such a spin flip process 
TBI . The toy model that was used in this study consisted of a particle with 



spin, having only a single translational degree of freedom, in the presence 
of a ID inhomogeneous magnetic field. We found that the trapped state of 

the particle decays with a lifetime given by ~ 1/ (^\/KuJmbj exp {2/K) where 

K = ujmb/^^prec, and where the result is valid for K <^ 1. Though the field 
that was used in this model did trap the particle, it was not realistic in 
the sense that it was not curl-free. Our next step was to study the case of a 
particle with spin, having two translational degree of freedom, in the presence 
of a physically more realistic trapping field that, in contradistinction to the 



toy model, is curl-free [|T9[. This model is reminiscent of a loffe-Pritchard 
trap [^, I, but without the axial translational degree of freedom. Here 
we found that the lifetime is given by ~ l/ci;„if,exp {2/K) which is similar to 
the result found in the ID case. In the present paper we describe an analysis 
of a loffe-Pritchard type trap which includes the axial translational degree 
of freedom. We neglect the effect of interactions between the particles in the 
trap, and we analyze the dynamics of a single particle inside the trap. Unlike 
in our previous papers, where we studied the case of a spin 1/2 particle, we 
treat here the case of a spin 1 particle, both as an example to show the validity 
of our approach for higher spins, and also because it is more relevant in view 
of the recent development in Bose-Einstein condensation experiments. 
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The structure of this paper is as follows: In Sec.(0) we start by defin- 
ing the system we study, together with useful parameters that will be used 
throughout this paper. Next, we carry out a classical analysis of the prob- 
lem in Sec.(^). Here, we find two stationary solutions for the particle inside 
the trap. One of them corresponds to a state whose spin is parallel to the 
direction of the magnetic field whereas the other one corresponds to a state 
whose spin is antiparallel to that direction. When considering the dynamical 
stability of these solutions, we find that only the antiparallel stationary solu- 
tion is stable, as expected from the discussion in Sec. fpT^ ) above. In Sec.(^ 
we reconsider the problem from a quantum-mechanical point of view for a 
spin-one particle. Here, we find states that refer to antiparallel (M = — 1, 
where M is the magnetic quantum number) and orthogonal (M = 0) orien- 
tations of the spin, the first of these being bounded while the second one is 
unbounded. We argue that the third possible situation, in which the spin is 
parallel (M = +1) to the direction of the field, has negligible coupling to the 
bound state, and therefore can be neglected. We show that the antiparallel 
and orthogonal states are coupled due to the inhomogeneity of the field, and 
we calculate the transition rate from the bound state to the unbounded state. 
Finally, in Sec.(^ we compare the results of the classical analysis with those 
of the quantum analysis and comment on their implications for practical 
magnetic traps. 



2 Description of the problem. 



We consider a particle of mass m, magnetic moment fi and intrinsic spin S 
(aligned with /i) moving in an inhomogeneous magnetic field B corresponding 
to traps of the loffe-Pritchard, 'clover-leaf and 'baseball' type 0, and given 
by 



B 



Bo + Ib"z' - \b" (x^ + y') 



(2) 



+ [B' - ^B"z I XX- 



-B' - ^B"z j yy. 



This field possesses a nonzero minimum of amplitude at the origin, which is 
the essential part of the trap. The Hamiltonian for this system is 

p2 

2m 



H 



/xB 



(3) 
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where P is the momentum of the particle. 

We define Up as the precessional frequency of the particle when it is at 
the origin = 0, y = 0, ^ = 0). Since at that point the magnetic field is 
B =Boz we find that 



IjlBq 



(4) 



Next, we define uJz and Ur as the small-amplitude axial and lateral vibrational 
frequencies of the particle when it is placed with antiparallel spin into the 
adiabatic potential given by 



ye// = /x|B|~/x5o(l 

from which we get 



BqB" 
2BI 



{B' f - '^BoB" , ^ ^ 



2BI 



UJr 



'liB" 



m 



/i [iB'f - iBoB"] 



mBn 



(5) 



In what follows we assume that {B') — \BqB" > such that Ur is real. We 
also define the ratios, 



B"S^ 
finiBQ 



(6) 



[B'f - \BoB"] 



ImuBq 



These parameters will be our 'measure of adiabaticity'. It is clear that as 
and Kj. become smaller and smaller, the adiabatic approximation becomes 
more and more accurate. Note that when the bias field Bq vanishes, both 
Kz and Kr become infinite, and the adiabatic approximation fails. We will 
show below that, under this condition, the system becomes unstable against 
spin flips, which is in agreement with our discussion at the beginning. This 
shows that the introduction of the bias fleld Bq, is essential to the operation 
of the trap with regard to spin-flips. 
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3 Classical analysis. 



3.1 The stationary solutions. 

We denote by n a unit vector in the direction of the spin (and the magnetic 
moment). Thus, the equations of motion for the center of mass of the particle 
are 

d'^x d 

'^^=%^''-^^ 

d 

and the evolution of its spin is determined by 

S— = fihxB. (8) 

(JjL 

The two equilibrium solutions to Eqs.(^ and (|) are 

h{t) = Tz (9) 

with 



x{t) = 
y{t) = 
z{t) = 

representing a motionless particle at the origin with its magnetic moment 
(and spin) pointing antiparallel (n(t) = — z) to the direction of the field at 
that point and a similar solution but with the magnetic moment pointing 
parallel to the direction of the field (n(t) = +z). 
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3.2 Stability of the solutions. 

To check the stabihty of these solutions we now add first-order perturbations. 
We set 

n(t) = Tz+e^{t)^+ey{t)y (10) 
x{t) = + 6x{t) 
y{t) = + 6y{t) 
z{t) = + 6z{t) 

(note that, to first order, the perturbation Sh = ex{t)'x.+ey{t)y is orthogonal 
to the vector n for the stationary solution no = =Fz, since n is a unit vector), 
substitute these into Eqs.(0) and (||), and retain only first-order terms. We 
find that the resulting equations for 6x{t), Sy{t), 6z (t), e^it) and ey{t) are 

m^ = ±^fiB"6x + fiB'e, (11) 



(f5y ,1 ^ 



S-nr = -f^Boe^ =F fiB'5x. 
at 

The motion of the 2;-coordinate is decoupled from the others. If B" > 0, it is 
stable only when the upper sign is taken, corresponding to a spin antiparallel 
to the direction of the field. It can be shown that when B" < then, even 
if the system is stable under axial vibrations (by choosing the lower sign), it 
cannot be stable as a whole. We therefore disregard the lower sign, and the 
equation for the z-coordinate for the rest of the derivation. 

The normal modes of the reduced system transform as the irreducible rep- 
resentations of the symmetry group. The 4-dimensional linear space spanned 
by the deviations {5x, Sy, ex, e^) from the stationary state carries the irre- 
ducible representations r+ with characters e"*''' and r_ with characters e'^^'^, 
and may thus be decomposed into the two 2-dimensional invariant subspaces 
transforming as r+ and r_, respectively. These subspaces are spanned by 
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the circular position coordinates and precessional spin coordinates 

T+ : {p+ = 6x + i6y,e- = ex - ity); (12) 
T_ : {p_ = Sx - i6y,e+ = + iey). (13) 

Thus, the normal modes consist of a circular motion in the {x,y) -plane cou- 
pled to a precession of the spin vector in the opposite sense. 

Indeed, after introducing the (p±, e^)-coordinates into Eqs.(|TT]), this set 
of four equations decomposes into one pair of equations for (p+,e_) and 
another pair for (p_,e+). We now look for oscillatory (stable) solutions of 
these equations and set 

(14) 

) , (15) 
) • (16) 

These equations have non-trivial solutions whenever the determinant of either 
of the two matrices vanishes. This yields the secular equations 

which determine the eigenfrequencies u of the various modes. Since the 
reduced system has three degrees of freedom, we expect to have three normal 
modes. Indeed, when is a solution of the first equation, then —u is a 
solution of the second equation. We define the mode frequencies in Eq.(p!7D 
to be positive (or, in the case of complex u, to have positive real part); 
the negative cu-values are needed to construct real solutions. Then, the r+- 
modes describe vibrational motions turning counter-clockwise coupled to spin 
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This yields the algebraic equations 

1 
2 



ipB' i {cuS + fiBo) J V e->o 



Y . I + mu^ fiB' \ f p-,0 



-ipB' I {uoS - pBo) J \ ^+.0 



precessions turning clockwise, i.e., opposite to the natural spin precession, 
and the r_-modes describe vibrational motions turning clockwise coupled 
to spin precessions turning counter-clockwise, i.e., in the same sense as the 
natural spin precession. 

Stability requires that all three solutions of, say, Eq.([T8|) be real. We 
note that at the edge of the stability region (and when Kr 7^ 0), two out 
of the three roots of Eq . (p^Sf) for uj become identical. In this case, the third 
order polynomial Eq.(|l8D takes the form P (u) = (cu — cJi)^ (tu — UJ2), which 
satisfy dP / du!\^^^_^ = 0. The edge of the stability region is then found by 
simultaneously solving the equations P (u) = and dP/du = 0. The result 
is given in the form of the parametric curve in the {K^,K^)-p\&ne 



which is shown in Fig.(|l|). Note that by eliminating t from the second equa- 
tion and substituting it in the first gives explicitly in terms of K"^. 

4 Quantum-mechanical analysis. 

4.1 The Hamiltonian and its diagonalized form. 

In this section we consider the problem of a neutral particle with spin one 
{S = h) in a 3D inhomogeneous magnetic field from a quantum- mechanical 
point of view. Unlike the classical analysis, in which the derivation was valid 
for any value of the adiabaticity parameters Kr and K^, we concentrate here 
on the behavior of the system when Kr and Kz are small. Note also that, 
quantum-mechanically, the magnetic moment /i and the spin S* of a particle 
are related by 



where 7 is the gyromagnetic ratio of the particle. Setting n = 7S' and S = h 
in Eqs.(|D gives 
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Now, it is convenient to transform to cylindrical coordinates (r, 0, z) by 
setting X = r cos (p, y = r sin 0. We denote by B the amplitude of B, by ^ its 
direction with respect to the z-axis and by cp the angle between the projection 
of B onto the (x,?/)-plane and the x-axis. Thus, Eq.(|D is rewritten as 

B = 5 [sin 6 cos (f'k+ sin 6 sin (fy + cos 6z] . (19) 

The approximate expressions for B, 6 and near the origin are given by 



B{r,^,z)c^Bo 1 



BoB" , , fiB')^-lBoB" 



2BI 



2BI 



4 2 2 

r ,zr ,zj 



(20) 



^ (r, 0, 2;) 



(y9 (r, 0, 2;) 



arctan 



arctan 



V'^J+5f\ B' 



r 



B, 



By 

B.. 



arctan — 



Bo 

y 



O {r^,rz,z^) 



X 



1 + 



B" 
'B' 



-0 + 0(2 sin (20)) 



Thus approximately, B and 9 depend only on r, whereas (p depends only 
linearly on 0. 

The time-independent Schrodinger equation for this system is 



— uB (sin 6 cos sin sin uDSy + cos Osz) 

2m 



E^{r, 0, z) 
(21) 



where s^, Sy and Sz are the spin one matrices, given by 
1 





-i 













i 


—i 












i 


j 




^ 








E is the eigenenergy, and $ is the three-components spinor 

$ = 



*+(r,0, 2;) 

<l>o(r,0,z) 
$_(r,0, z) 



(22) 



In matrix form Eq.(PTD becomes 

(Hk + Hm) I $o(r,0,2) 
$_(r,0,2;) 



<!>+{r,^,z) 
^4r,(j),z) 



(23) 
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where and Hm, given by 



2m 



I 



(24) 



Hm = — /i-B 



cos 6* 

isin^e*'^ 



1 







sm ( 



-up 



72 



sm I 



V2 
— cos 6 



J 



are the kinetic part and the magnetic part of the Hamiltonian if, respectively. 

In order to diagonahze the magnetic part of the Hamiltonian, we make 
a local passive transformation of coordinates on the wavefunction such that 
the spinor is expressed in a new coordinate system whose z axis coincides 
with the direction of the magnetic field at the point (r, 0, z) . We denote by 
R (r, 0, z) the required transformation and set \l/ = Thus, ^ represent 
the same direction of the spin as before the transformation but using the 
new coordinate system. The Hamiltonian in this newly defined system is 
given by RHR~^. We represent the rotation matrix R in terms of the three 
Euler angles: First, we perform a rotation through an angle ip around the 
z axis. Second, we make a rotation through an angle 9 around the new 
position of the y axis. At the end of this process the new z axis coincide 
with the direction of the magnetic field. Now the value of the last Euler 
angle, which is a rotation around the new z axis, has no effect on this axis. 
For simplicity we choose this angle to be 0. Thus, the representation of the 



complete transformation for spin-one particle is given by ||22 



R = exp [iOsy] exp [iipsz] , 

while its inverse is given by 

R^^ = exp [—i(psz] exp [—iOsy] . 

It is easily verified that the transformation indeed diagonalizes the magnetic 
part of the Hamiltonian as 



RHmR 



-fiBsz- 
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As for the kinetic part we show at the Appendix that 



RHj^R ^ = 



2m 



—iV'^(f {cos 6sz — sin^s^,.) — jVy^l^ (cos^^s^ — sin^^s^)^ 
—2i (cos 6*52 — sin 6*53;) Vv? ■ {—iV9sy + V) 



iv^eSy - 1 v^r si - 2isyVe ■ v + 



(25) 



Since we are interested in the behavior near the origin, we substitute the 
approximate expressions Eqs.(pOD into Eq. (p5D , replace cos^ by 1 and sin^ 
by (since 6 changes very slowly as compared to the extent over which fiB 
changes significantly), and neglect the terms that are proportional to 
and |V^|^ (being of higher order with respect to ■ V ). This gives 



RHkR ^ 



2m 



2is, 



B' d 
Bq dr 



1 

2«^s, 



d 



^s. 



Thus, the Hamiltonian of the system in the rotated frame may be written 
approximately as 



H ^H, 



diag 



Hi, 



where 



Hdiag 



H 



2m 
h^B' 



int 



d 



mBg dr 



2i ^ d 



inti 



^S, 



- l^Bsz 



(26) 



(27) 



The first part of the Hamiltonian Hdiag is diagonal with respect to the spin 
degrees of freedom. It contains the kinetic part ~ V^, a term whose form is 
—fiBsz which is identified as the adiabatic effective potential, and the terms 
~ 1 / r"^ , ir~'^ s zd / d(j) which appear due to the rotation. The second part of 
the Hamiltonian Hint contains only non-diagonal components. Generally, 
Hint should contain terms which couple a spin state M to the two nearest 
spin states M ± 1 and to the two next-to-nearest spin states M ± 2 (see 
the Appendix). In the limit where Kz and are small, we see that the 
coupling of the state with spin projection value M to the states M ± 2 is 
negligible compared to its coupling to the Mil states. We proceed to find 
the eigenstates of Hdiag- 
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4.2 Stationary states of H^iag' 

Since Hdiag is diagonal, the three spin states of the wavefunction are decou- 
pled. We seek a solution for the spin-down (M = — 1) state 

/ 

^_ = \ ; E = E_, (28) 

\ 'ip_{r,(f),z) 



and for the M = state 

^0 





ilJo{r,(j),z) 




E = En 



(29) 



We do not consider the spin-up (M = +1) state, since its coupling to the 
trapped spin-down state is negligible, as explained above. 

The equation for the non-vanishing component of the spin-down state 
reads 



2m 



(30) 



whereas the equation for the non- vanishing component of the spin-zero state 



IS 



2m 



(31) 



The solutions of these equations is outlined in the next two subsections. 



4.2.1 Stationary spin-down (M = —1) states. 

Eq.(^) represents a particle in a cylindrically symmetric attractive 3D po- 
tential. If the extent of the wave function is small enough we can expand 
B in Eq.(pU[) to second order in r and z as given by Eq. ([2(]|) , and apply 
the well-known solution of the harmonic oscillator in 3D. Under this 
approximation, Eq.(^0|) becomes 




r 




(32) 
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The z-coordinate decouples and we assume that it is in the ground-state. We 
thus seek a solution whose form is 



V^_(r,^,0) = /(r)e^^^f 



1/4 



exp 



mw^z 



2h 



(33) 



with u an integer. The equation satisfied by / (r) is then 



2m 



d'f Idf f 2' 



+ 



1 1 
mujtr 



f= {E_-fiBo--nLuJf. 



(34) 



This is an eigenvalue problem for /. The smallest eigenvalue for this equation 
is obtained by setting 



1/ = 1 

for which the eigenfunction / is 

/ (r) = De'"^ exp 



r 2 



2h 



Thus, under the harmonic-oscillator approximation, the normalized down- 
part of the spin-down state is given by 



e ^ exp 



mujr 



^'2 



2h 



exp 



niLU 



^ ^2 



2h 



(35) 



Note that the extent of this wave function over which it changes appreciably 
is given by 



/ B 

Az ~ x/kI\ — ; Ar 




(36) 



whereas the extent over which changes significantly (see Eq.(pOD) is 

^r,B - M, (37) 



Thus, the ratio between these two length scales is 



Az 

AVf^B 



Ar 

AVf^B 



'Kr. 



(38) 
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We therefore conclude that when and are small enough, the harmonic 
approximation is justified. 

The wave function ?/^_, given by Eg. (|35|) , then represents the lowest possi- 
ble hound state for this system. This state corresponds to a trapped particle. 
The energy of this state is 



= /i5o + -jf^z + thJr = /i5o [1 + Kz + 2Kr) ~ /i-Bo 



while its full spinor representation is 



v 








e "^ exp 



TT 



h 



2h 



exp 



muj 



2h 



•2 2 



(39) 



(40) 



4.2.2 Stationary (M = 0) states. 

Eq.(^) describes a free particle. It corresponds to an unbounded state rep- 
resenting an untrapped particle. In this case there is a continuum of states, 
each with its own energy. As we are interested in non-radiative decay, we 
focus on finding a solution with an energy which is equal to the energy found 
for the trapped state, that is 

Eo = E_^ liBo. (41) 
We seek a solution in the form 

ipo^r, (p) = g{r) exp [ik^z + , 

where (3 is an integer. Substituting this, together with Eq. f^I]) into Eq. (^T|) 
gives 

_dr'^ 

where 



The non-singular solution 



.11. (*?-^ 

r dr \ 



2iJimBQ 



for g is 

g (r) = Jp (Kr) , 
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where Jp (x) is the Bessel function of the first kind of order j3. For what 
follows, it is convenient to introduce an angle 7 such that 



kr = kf) sin 7 
kz = ko cos 7. 



/ 2/imSo 

with < 7 < vr. 

We note that Hint does not operate on the coordinate. Hence, in order 
to have a non- vanishing matrix element between the zero-state and the down- 
state, they must have the same 0-dependence. Thus, /3 = z/ = 1, and as a 
result, the state with angle 7 is given by 



'4'Q{r, 0, z) = C^Ji {kor sin 7) exp [i (0 + k^z cos 7)] . (42) 



with 





\E'o = I CjJi {kor sin 7) exp [i (0 + k^z cos 7)] | , (43) 



where C-y is the normalization constant which is chosen to be real, and de- 
pends on 7. To evaluate Cj we temporarily introduce boundary conditions 
under which the wavefunction \E'q vanishes at r = R, and satisfies periodic 
boundary conditions along z with period Z. Thus, normalization of gives 



/Z/2 p2TT pR -t 

dz dct) rdr\ml\^ = C^^Z2TX-R^[J2{hR?,m-i)f 
z/2 Jo Jo 2 



such that 



^R\J2{koRsm-f)\ 



where we have used 



r[Mkr)frdr = l-R'[J2ikR)f. 
Jo ^ 
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In the asymptotic region kR > 1, the function J2{kR) takes the values 



±A/2/(7rA;i?) at the zeros of Ji{kR). Thus, 

1 a/^o sin 7 



/Z^R\J2ikoRsm-f)\ V2ZR 
and hence 

^ (44) 
2ZR ^ ^ 

4.3 The transition rate. 

We calculate the transition rate from the bound state given by Eg . (^OD to 
the unbounded state Eg. (|i3|) , according to Fermi's golden rule Thus, 
the infinitesimal decay time from the trapped state to the untrapped state 
defined by 7 is given by 

^(^) =YI^'|' -^7(^0)^7, (45) 

where p^{E)d'~f is the density dN/dE of states ipl with an angle between 
7 and 7 + (i7 and energy between Eq and Eq + dE, and if/ is the matrix 
element of Hint between the bound state and the unbounded state defined 
by 7 and Eq. To find p^{Eq) we note that the final state is defined by the 
two guantized /c- vectors kr = k^sm'^ and k^ = /C0COS7. The possible kz 
values are egually-spaced with lattice constant dk^ = 2n/Z. Since the Bessel 
function Jiikrv) is very close to its asymptotic behavior at large arguments 



Ji{krR) ~ \ cos IkrR - ^ 

V nkrR \ 4 

the kr are also very much egually-spaced (even when kr is small, it is still a 
good approximation) with lattice constant dkr — n/R. Thus, in the {kr, kz)- 
space, the allowed k vectors form a regular lattice, and the number of states 
dN in the volume element k^dkod'y is given by 

kod'jdko 

dN 



\RJ V Z 
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With dE = h'^kodko/m, this gives the density of states 



dE 27r2^2 



This, together with Eq.(0) yields 



mko sin 7 



Evaluation of gives 

00 2lT 

poo 

H7 = 



dz jrdr j dcj)^]^ Hint^ . 



27r 







X / drr'^Ji [k^r sin 7) exp 



. WB' (-i) 



2/i 



-r 



muJr 
dz exp 



(46) 



^ z^ — ikoz cos 'y 



2h 



and hence 



If/' = 47r^- 

7ry 



X exp 



\V2mBo) 
hkl cos^ 7 /cq sin^ 7 



h y nh 



/mciJrX^ ^2 ^tt/i /A;o^^sin7Y 



(47) 



Substituting Eqs. (^Tj) and ( ^^D into Eq.(^5D and integrating over 7 from to 
71 gives 



(48) 



with 



/•TT 

/ (a, h) = I d'y sin^ 7 exp [—a sin^ 7 — 6 cos^ 7] 
Jo 



(49) 
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where we have substituted our previous definitions for Up, oOr and 00 z- The 
integral / (a, b) may be expressed in terms of the simpler integral 



/•TT 

Iq (a, b) = d'y sin 7 exp [—a sin^ 7 — 6 cos^ 7] 
Jo 

= 2 exp (—a) / exp [(a — 6) t^] dt 
Jo 



(50) 



by 



/ (a, 6) 



da 



(51) 



In the isotropic case where Ur = ujz = the integral in Eq. (^8D can be 
evaluated analytically with the result that 



T 



Sv 271 y/ujpUJi exp 



In the extreme cases Ur ^ uj^ and Uz S> cc'r we obtain from the asymptotic 
behavior of the error function of real and imaginary argument 



Jo (a, b) 



exp (—a) ] b ^ a 



1 

- exp (—6) ; b <^ a 
a 



(52) 



and hence 



/ (a, 6) 



^/q (q, fe) 
da 



— exp (—a) ; b:^ a 
■ exp (—6) ; b <^ a 



(53) 



Substituting Eqs. (M) and (pSQ into Eq.(|g) gives 

2cj^ 



1 



T 



47ru;r exp 
8 ^/2tt y/cOpCOi exp 



2U!„ 





2ujp 








UJz _ 



for tUp ^ UJz 



; for cijp ^ ciJr = = 
; for cijp ^ cj^ ^ 



(54) 



with the conclusion that the transition rate is dominated by the largest of 
the two frequencies Uz and Uj.. 
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5 Discussion. 



The problem we have studied has three important time scales: The shortest 
time scale is Tprec, which is the time required for one precession of the spin 
around the axis of the local magnetic field. The intermediate time scale is 
given by Tr^z = Tprec/ K^^z, which are the times required to complete one cycle 
of the center of mass around the center of the trap in the lateral and axial 
directions, respectively. These two time scales appear both in the classical 
and the quantum-mechanical analysis. The longest time scale (provided that 
Kj. and Kz are small) T^sci which is not present in the classical problem, is 
the time it takes for the particle to escape from the trap. 

Whereas the classical analysis yields an upper bound for Kz and for 
trapping to occur, no such sharp bound exists in the quantum-mechanical 
analysis. Nevertheless it is interesting to compare the classical bound with 
the values of Kz and K^ for which the exponent in the expression for the 
quantum-mechanical lifetime becomes equal to 1: According to Fig.(|l]), we 
find that ir^,max = l/v^ = 0.707 when Kr = 0, and /sT^^max = ^4/27 = 0.385 
when Kz = 0. From Eq.(|5^) on the other hand, we conclude that Kz,cr = 0.5 
when Kr = 0, and Kr^^r = 0.5 when Kz = 0. Thus, the quantum-mechanical 
condition for trapping to occur is roughly the same as the classical condition. 
These results however, should be taken with caution since our quantum- 
mechanical analysis is valid only for small values of Kr and Kz- 

Though our derivation was for the case of a spin-one particle, it is clear 
that it can be extended to particles with higher spin, and also to half-integer 
spin particles. In view of the results obtained by our recent study of spin 
half particles in ID field |jl8| and 2D field |]19|, we believe that the expression 
for the lifetime in these cases is similar to the result which is obtained in the 
present paper. 

As an example, we apply our results to the case of a spin 1 atom that is 
trapped in a field with Bq = 100 Oe and Bq/B' ~ ^Bq/B" ~ 10cm. These 
parameters correspond to typical traps used in Bose-Einstein condensation 
experiments ||2^, |3^ . The results, being correct to within an order of 

magnitude, are outlined in Table |I]. We note that in both cases the values 
of Kz and Kr are much smaller than 1. Also, the calculated lifetime of the 
particle in the trap is extremely large, suggesting that the particle is tightly 
trapped in this field. 

In this study we have been interested in the ground-state trapped state. 
In the case of a particle with spin 1/2 or spin 1, this is the only one trapped 
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spin state. However, when particles with higher spin are considered there are 
more than one trapped states. A natural question in connection with these is 
what is the lifetime of these trapped states. Another interesting issue is the 
lifetime of an excited state in a given trapped spin state. Our preliminary 
results show that some of these excited states may have a short lifetime, 
being algebraically dependent on cUp/cUr and cUp/cUz rather than exponentially 
dependent. This question is still under study. 
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A Transformation of V^. 



The transformation of is given by 

RV^R-^ = exp [iOsy] Q exp [-iOsy] , (55) 

where 

Q = exp [i(pSz\ exp [— ic/^s^] . (56) 

Evaluating Q first gives 

(exp [-i(psz] A) = V ■ V (exp [-iips^] A) 

= V ■ [V (exp [-iiph]) A + exp [-i(ps^] VA] 
= V ■ [—i exp [—iLpSz] VipSzA + exp [— i(/9S2] VA] 
— V • [—i exp [—iipsz] Vips^A + exp [—iipSz] VA] 

but 

V • [—i exp [— ic/js^] V(pSzA] 

— Vip • V [—i exp [—icpsz] SzA] — i exp [— if^s^] V'^cpszA 

— V(fi ■ [—i exp [— ic/js^] SzVA — exp [—i(fiSz] V<^s^>l] — i exp [—iipsz] V'^ipSzA 

and 

V ■ [exp [—iipsz] VA] 

— —iVip exp [—iipsz] Sz ■ VA + exp [—iipSz] V'^A 

hence 

(exp [-iifSz] A) 

— V</7 • exp [— i(/?S2;] A — exp [— V<^s^>l] — i exp [— ic/js^] V^f/Ps^A 

— iV</?exp [— • VA + exp [— 

= — 2iV</?exp [— ic/ps^;] • VA — exp [— if/ps;^] |V</?|^ s^A 

— i exp [— V'^(fiSzA + exp [— ic/js^] V^>1 

thus 

exp [ifSz] (exp [—iipsz] A) 

= -2iszVip ■ VA - \V(pf s\A - iV'^ipszA + V^A 
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or in an operatorial form 



Q = exp [iifSz] exp [—iifSz] = —2iSzVif ■ V — \V(ff — iV'^ifs^ + V^. 



(57) 



Substituting in Eq.(|55|) each of the four terms in Eq.(|57D we find 



exp [iOsy] exp [-iOsy] = -iV'^Osy - | V^|^ si - 2iSyVe ■ V + V^ (58) 



exp [iOsy] (— zV^v^) s^exp [—iOsy] = —iV'^ip {cosOsz — sin6'Si,.) , (59) 



exp [iOsy] si exp [—iOsy] = (cos 9sz — sin Os^ 



(60) 



cos^ 6sl — sin 6 cos 6 is^s^ + s^Sz) + sin^ 6s 



2 /3s2 



exp [iOsy] SzV exp [—iOsy] = exp [iOsy] s^ exp [—i9sy] {—iVOsy + V) (61) 

= (cos6's2 — sin6'sK) {—iVOsy + V) 
= —i cos BVOszSy + i sin OVOs^Sy + cos V — sin 6's2; V, 

Substituting Eqs.([58|) to (|61|) into Eq.(|55D gives 



2 n-1 



— iV^V9 (cos 6^52 — sin6's^) — \ {cos6sz — sin6'Sa 



—2i (cos 9sz — sin 6^5^) Vip ■ {—iVOsy + V) 

-iV^esy - \vef si - 2isyVe ■ v + 

Note that the transformed is composed of terms containing s" with n = 
0, 1 or 2. This is a consequence of the fact that the original operator 
is a second order differential operator. Thus, a spin state \1/a/ for which 
■Sz^&M = M'ifM is coupled, in first order, only to the states ^^Mii and \&m±2- 
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Table 1: Typical time scales for a spin 1 atom trapped with a field Bq = 100 
Oe and Bq/B' ~ y^Bp/B" ^ 10cm. 





Spin 1 atom 


m gr 


~ 10-22 


fj, emu 


~ 10-20 


K., Kr 


~ 10-8 


uj~^ sec 


~ 10-9 


uj~^, uj~^ sec 


~ 10-1 


T^sc sec 


~ io(^°^) 
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